Critical properties of loop percolation models with optimization constraints 
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We study loop percolation models in two and in three space dimensions, in which configurations 
of occupied bonds are forced to form closed loop. We show that the uncorrelated occupation of 
elementary plaquettes of the square and the simple cubic lattice by elementary loops leads to a 
percolation transition that is in the same universality class as the conventional bond percolation. 
In contrast to this an optimization constraint for the loop configurations, which then have to 
minimize a particular generic energy function, leads to a percolation transition that constitutes a 
new universality class, for which we report the critical exponents. Implication for the physics of 
solid-on-solid and vortex glass models are discussed. 
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I. INTRODUCTION 



The percolation of loops (or closed strings) appears 
naturally in the context of liquid helium |l], ^] , early uni- 
verse j| || and high-temperature superconductors ||, [Hs) . 
where loops represent world- lines, cosmic strings and vor- 
tex loops, respectively. In analogy to the characteristic 
size of a cluster in conventional site or bond percolation 
H , the typical diameter £ of the loops diverges when ap- 
proaching a critical point, the loop percolation transition. 
This transition shows power-law behavior at the critical 
point, which is described by a set of critical exponents 
that constitute a universality class. 

In this paper we study these percolation transitions 
numerically. We show that the loop percolation (LP) 
model, in which each elementary plaquette of a square 
lattice (in 2d) or on a simple cubic lattice (in 3d) is occu- 
pied with a probability p with an elementary loop, is in 
the same universality as the conventional d-dimensional 
bond percolation. More importantly we also show that 
in contrast to the LP model an optimization constraint 
for the loop configurations, which then have to minimize 
a particular generic energy function, leads to a percola- 
tion transition that constitutes a new universality class, 
for which we report the critical exponents. 

These loop Hamiltonian (LH) models are relevant for 
the ground state properties of disordered solid-on-solid 
models Jll], [I2J [l3| and vortex glasses jl4| . One particu- 
lar example of was recently studied by us in the context 
of a 3d vortex glass model for amorphous high-T c super- 
conductors in the strong screening limit |TH] , which shows 
an unconventional percolation transition of vortex loops 
in the ground state as a function of the disorder strength 
a. 

The paper is organized as follows: In section [0], we in- 
troduce the models and the definition of loops (clusters). 
In section |l| we locate the percolation transition in each 
model and calculate the critical exponents using finite- 
size scaling (FSS). Section [V concludes the paper with 



II. MODELS 

Consider a <i-dimensional hyper-cubic lattice - i.e. a 
square lattice in 2d or a simple cubic lattice in 3d - of 
linear size L (L = 7 in the 2d example in Figjl]) with 
free boundary conditions (f.b.c). In the loop percola- 
tion (LP) model the elementary plaquettes of the lattice 
are occupied with elementary loops with a probability 
P\. (There are (L — I) 2 elementary plaquettes in 2d and 
3L(L — l) 2 elementary plaquettes in 3d.) An elementary 
loop consists of the four bonds belonging to an elemen- 
tary plaquette plus a randomly chosen direction: cither 
clockwise or counterclockwise both with probability 1/2 
(see Fig. |l|). When two adjacent plaquettes are occupied 
by elementary loops of same orientation we cancel the 
occupation of the common bond as indicated in Fig. [l| 
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a summary and a discussion. 



FIG. 1: Configuration n of the loop percolation model on a 
2d square lattice with system size L = 7. 

We can identify the resulting (directed) bond configu- 
ration of the LP model with a flow n = {ni,ri2, ...,tim}j 
where rii is an integer and M is the number of bonds 
in the lattice. We say that rii = if bond i is not oc- 
cupied, rii = ±1 if it is singly occupied in the positive 
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FIG. 2: Typical loop configurations of the LP model around the critical point pi c » 0.41 in 2d (top) for L = 50 and pi c » 0.08 
in 3d (bottom) for L = 8. In 2d the different loops are marked by different grey scales (colors), whereas in 3d all loops are 
black except for the percolating loop, which is marked by light grey (red). 



(negative) direction (positive and negative are defined by 
the introduction of an appropriate coordinate system), 
rii = ±2 if it is doubly occupied etc. Thus, an elemen- 
tary loop (e.g. in the xy plane) can be represented by 
n XtV = n x +x, y = —n x +i, y +i = -n x , y+1 = 1 if oriented 
counterclockwise, as shown in Fig. [jj. The complete flow 
n then can simply be thought of as the sum of all ele- 
mentary loops. Obviously, in this sum the flow variables 
on the common bonds of adjacent elementary loops can- 
cel arithmetically. Moreover, the construction of this flow 
via addition of elementary loops implies that on each site 
of the lattice the number of ingoing arrows balances the 
number of outgoing arrows (see Fig. |l|): one says that 
this flow is divergence free 

V-n l = 0. (1) 

(The lattice-divergence operator is defined on each lattice 
site and sums simply all 2 • d flow variables of the bonds 
connected to it). 

In analogy to conventional bond percolation we are 
now going to define the clusters of a configuration of the 
LP model: Two occupied bonds belong to the same clus- 
ter if they have one site in common. Thus, all bonds of 
the cluster can be connected via a directed path along oc- 
cupied bonds belonging to this cluster, which is analogous 
to conventional bond percolation - up to the attribute di- 
rected. This is actually a slightly non-trivial observation 
- see Fig. |l| to exemplify this statement - because an oc- 
cupied directed bond cannot be traversed in the opposite 



direction, but is a direct consequence of the fact that a 
cluster is also a sum of elementary loops. 

The mass m of a cluster is the number of occupied 
bonds, where a bond i with a flow n\ counts \rii\ times, 
i.e. has mass rrii = \rii\. A percolating cluster is a cluster 
spanning the entire system (in at least one of the d direc- 
tions). This implies that the cluster contains a directed 
path along its occupied bonds from one side of the lattice 
to the opposite. In the following we refer to the cluster 
just defined as a loop. 

The representation of the configuration of directed oc- 
cupied bonds as a flow n is now used to define the LP 
model with an optimization constraint: In contrast to the 
stochastic uncorrelated occupation of elementary loops in 
the LP model above, we now consider a different type of 
occupation of loops, which results from the minimization 
of an energy function for a loop configuration n 

M 

H = H{n} = J2Mn l ) , (2) 

i=l 

where the sum is over all bonds i on a d-dimensional 
hyper-cubic L d lattice with periodic boundary conditions 
(p.b.c). Herewith we assume the energy function to be 
composed of solely of local terms fi(n) with fi(n) > 
and convex (i.e. f"{n) > 0) for all bonds i and flow (or 
occupation) values n. Such energy functions are relevant 
in the context of disordered solid-on-solid models (in 2d) 
p"T| , O, |l3| and vortex glasses (in 3d) Il4j since they 
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FIG. 3: Typical loop configurations of the loop Hamiltonian (LH) model in the ground state (T = 0) around the critical point 
o c ~ 0.46 in 2d (top) for L — 50 and a c ~ 0.31 in 3d (bottom) for L = 16. In 2d and 3d all loops are black except for the 
percolating loop, which is marked by light grey (red). 



determine their ground states. One has to keep in mind, 
that the loop condition (1) has to be fulfilled — i.e. the 
optimization task consist in finding the minimum of (2) 
under the constraint (1)! 

For fi(n) — f(n) independent of the bond index the 
minimum is trivial: n — 0, i.e. no bond is occupied. Only 
if the minima of local cost functions vary from bind in- 
dex to bind index in a non-trivial way one can expect 
a non-trivial loop configuration. We assume a random 
distribution of these minima (at values &,) and restrict 
ourselves to a quadratic form of fa (n) around these min- 
ima: fi(rii) = (rii — bi) 2 , which means that we study the 
loop configurations (i.e. occupied bond configurations 
that fulfill (1)) that minimize 



H = ]T(n t - hf . (3) 



The random variables hi are uniformly distributed hi € 
[—2a, 2a] at a fixed disorder strength a G [0,1]. Here, 
the probability p to occupy a bond depends on the dis- 
order strength a. We refer to the LP model with an 
optimization constraint of Hamiltonian (||) as the loop 
Hamiltonian model. For the minimization of a free en- 
ergy (optimization) we restrict to the calculation of the 
ground state (T = 0) configuration n, which is a min- 
imum cost flow problem that can be solved exactly in 
polynomial time with appropriate algorithms M . 



III. RESULTS 

We use a depth-first search algorithm known from 
combinatorial optimization Q| to identify the connected 
loops. The number of realizations we used to get statisti- 
cally averaged data varied from 500 for the largest system 
size to 20000 for the smallest system size. In the follow- 
ing the error bars of our data in the figures are smaller 
than the symbol size and are therefore omitted. 

Fig. U and || depict three typical loop configurations 
of the LP and LH model around the critical threshold, 
respectively, which clearly indicates a percolation phase 
transition for both models. 



A. Loop percolation model (LP) 

First, we study the LP model and consider the prob- 
ability pp erc that a loop percolates the system. Since 
we assume to have only one typical length scale, which 
diverges at the critical point like £ ~ \p% — p\ c \~ u , in a 
finite system pp £rc is expected to scale like 

P* a -°(L)~P[(p 1 - P i B )L 1 ' v ]. (4) 

Thus, pP erc ( K L) is independent of L at p\ c and the 
data curves should intersect for different system sizes 
L. From our raw data in the inset of Fig.^ (left) we 
locate the critical point p\ c = 0.4070 ± 0.0005 in 2d and 
p lc = 0.0793 ± 0.0004 in 3d, respectively We plot the 
scaling assumption (||) in Fig. || (left) and estimate the 
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FIG. 4: Finite-size scaling (FSS) for the LP model in 2d (top) and 3d (bottom): Plot of the percolation probability pP erc 
(left) and of the probability for a bond belonging to a percolating loop (middle). The inset shows the raw data. (Right) 
Plot of the average number n m of loops of mass m per lattice bond at pi c = 0.407 in 2d and at pi c = 0.0793 in 3d, respectively. 



inverse correlation length exponent \ jv — 0.75 ± 0.03 in 
2d and 1/v = 1.143 ± 0.090 in 3d from the best data 
collapse at fixed p ic . 

To get a second critical exponent we consider the prob- 
ability Pqo that a bond belongs to the percolating loop, 
i.e. the order parameter, which is expected to obey 



P oc (L)~L-P/»P[(p 1 -p lc )L 1 ^]. 



(5) 



Fig. U (middle) shows the raw data (inset) and the plot 
of the scaling law (|) with /3/v = 0.104 ± 0.020 in 2d 
and j3/v = 0.49 ± 0.02 in 3d such as to achieve the best 
data collapse. From the v above we determine (3 shown 
in Table § 

At the critical point p c \ the average number n m of 
finite loops of mass m per lattice bond scales like 



t(L,pi =pi c ) 



(0) 



where r is the Fisher exponent |6j, |8|. Since we assume 
the usual scaling relations of conventional percolation to 
be valid || , we also expect a combination of them to be 
valid, i.e. the hyperscaling relation 



2 - PHdv) 
1 - 0/(dv) ■ 



(7) 



where d is the spatial dimension. From the fit of the 
data of n m (L) at p cl in Fig. || (right) we get r depicted 
in Table || This is consistent with the value from putting 



the above v and (3 into Eq. (|7|), i.e. r = 2.05 ± 0.05 in 
2d and r = 2.20 ± 0.06 in 3d. 

To determine the critical probability p c that a bond is 
occupied we calculate 



p c (L) = m n m (L;pi = p lc ), 



(8) 



m— 4 



where mj, is the largest finite loop. We plot p c (L) versus 
1 / L as depicted in Fig. || and extract p c from the limit 
L — > oo. 

In addition to the results presented so far, we found 
that the mean number N perc of percolating loops per 
sample in the finite system can be described by a smeared 
step function with an upper boundary N perc = 1 as 
known from conventional percolation Q . When we define 
the mass wii of an occupied bond i to be rrii — 1 even for 
\rii\ > 1, the critical scaling behavior remains unchanged 
and the critical probability becomes p c = 0.565 ± 0.005 
in 2d and p c — 0.266 ± 0.005 in 3d, respectively. We 
also studied the case, where the algorithm detects the 
loops along oriented paths, and found the same results, 
as expected from what we said above. 



B. Loop hamiltonian model (LH) 

For the loop Hamiltonian (||) we perform analogous 
data analysis. From the intersection of the L-dependent 
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FIG. 5: FSS for the loop Hamiltonian (LH) model in 2d (top) and 3d (bottom). Plot of the percolation probability p perc 
(left) and of the probability P^ for a bond belonging to a percolating loop (middle). The inset shows the raw data. (Right) 
Plot of the average number n m of loops of mass m per lattice bond at cr c = 0.458 in 2d and at a c = 0.3129 in 3d, respectively. 



curves we locate 

the critical disorder strength at er c = 0.458 ± 0.001 in 2d 
and cr c = 0.3129 ± 0.0005 in 3d, respectively. From the 
finite-size scaling behavior of the percolation probability 
pperc ( similar t0 Eq ^ (|)) we get \ j v= 0.30 ±0.05 in 2d 

and 1/v = 0.95 ± 0.05 in 3d. The resulting exponent v 
are given in Table B. In 2d we find a value v = 3.3 ± 0.3, 
which is rather large. 

Fig. |B| (middle) shows the plot of the raw data of P^ 
(inset) and its scaling law similar to Eq. (^) with j3/u — 
0.55 ± 0.05 in 2d and (3/v = 1.30 ± 0.05 in 3d. From the 
v above we determine (3 in Table Q. 

In Fig. [5] (right) we plot the loop distribution n m (L) 
vs. the mass m at the critical point a c and determine r 
by power law fit, which gives the values shown in Table Q. 
From v and (3 above we get via the hyperscaling relation 
(0) the Fisher exponent r = 2.38 ± 0.17 in 2d and r = 
2.76 ± 0.26 in 3d, which are consistent with the values 
from the power-law fit within the error bars. 

In Fig. we plot the mean number N perc of percolat- 
ing loops per sample. The different curves of N perc (L) 
intersect at the same critical point a c found for pp erc 
above. Similar to the scaling law of N p ^v {L) m conven- 
tional percolation || we expect N perc (L) to obey 



N perc (L) ~ N[(a - a c ) L 1 ^} 



(9) 



and estimate the same v as above from the best data 
collapse, see Fig. [?]. 
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FIG. 6: Plot of the critical probability p c (L) vs. inverse sys- 
tem size 1/L in 2d and in 3d for the conventional percolation 
model, the LP model and the loop Hamiltonian (LH) model. 



The mean number N perc of percolating loops per sam- 
ple can become larger than one slightly above the critical 
point a c in contrast to conventional percolation, where 
only one percolating cluster exists for p > p c . The ap- 
pearance of several percolating loops can possibly be re- 
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for the loop Hamiltonian (LH) model in 2d (top) and 3d (bot- 
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FIG. 8: Schematic mapping of the loops in the LP model (left) 
onto sites (filled circles) in the conventional site percolation 
model with next nearest neighbors (right; see text). 



lated to the fact that the loop density in the LH model at 
a c is much smaller than in the LP model, as can be seen 
from comparison of typical loop configurations (c.f. Fig. 
U and [|). Moreover, the maximum of N perc (L) seems to 
increase with increasing L. From our data we could not 
determine the behavior of the maximum of N perc (L) in 
the thermodynamic limit L — > oo, in particular whether 
it converges to a constant or diverges. We also checked a 
L-dependent power law behavior of N perc (L) similar to 
Eq. (H) with a new critical exponent x (instead of 0) and 
found x/v = 0.01 ± 0.01 in 2d and x/v = 0.15 ± 0.02 in 
3d, i.e. x close to zero. 

Finally, we calculate the probability p c (L) that a bond 
is occupied (analogous to Eq. ||) for different L at the 
critical point a c as shown in Fig. |^ and extract p c as the 
average value, see Table |. 

In all considered ground state configurations of the 
loop Hamiltonian (^|) we found a bond to be empty or 
singly occupied only. Due to this observation we also 
investigated - beside the study presented this paper - a 
modified LP model, in which a plaquettc is allowed to be 
occupied with pi if and only if the amount of the result- 
ing flow n is |rti| < 1. Here, the algorithm checked each 
plaquette to be occupied or not in positional order. Note 
that this occupation process depends on the algorithmic 
order of occupying plaquettes in the system. Again, we 
found the same critical exponents as known from con- 
ventional percolation, but a different critical probability 

Pic- 

We also studied the 2c? LH model with a different prob- 
ability distribution function P(bi), where bi is given by 
a sum of two uniformly distributed random numbers out 
of [0, 2a]. This corresponds to the solid-on-solid (SOS) 
model on a disordered substrate, which has been stud- 
ied |ll|, 12 only at a = 1/2 yet. For this probability 
distribution function we get the same critical exponents 
as found above, but with a different critical point at 
a c = 0.395 ± 0.005, i.e. p c = 0.34 ± 0.02. This implies 
that our study is relevant to describe a disorder-driven 
flat-to-superrough phase transition, not studied in liter- 
ature yet. 

Closely related to the SOS model is the 2d model of 
a random elastic medium with contour loops, for which 
Zeng et al. found [|ll| the geometrical exponents f3/v — 
d-d f = 0.54 ± 0.01 and r = 2.32 ± 0.01 at a = 1/2. 
These exponents agree with the critical j3/v — 0.55 ±0.05 
and r = 2.38 ± 0.17 we found here at ov w 0.458. 



IV. SUMMARY 

We studied two loop percolation models, numerically: 
in the loop percolation (LP) model the loop configura- 
tion resulted from an uncorrelated unbiased random oc- 
cupation of elementary directed plaquettes, while in the 
loop Hamiltonian (LH) model the loop configuration ap- 
peared according to the Boltzmann weight of a particular 
microscopic model at T = 0, i.e. from an optimization 
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0.458 ±0.001 
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0.70 ± 0.02 


0.59 ±0.02 


0.64 ±0.02 


0.67 ±0.02 
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4/3 = 1.3 


1.33 ±0.05 


3.33 ±0.30 


1.33 ±0.05 


1.33 ±0.04 




5/36 = 0.138 


0.139 ±0.030 


1.80 ±0.35 


0.138 ±0.027 


0.139 ±0.007 


T 






2.45 ± 0.05 


2.05 ±0.10 




3d p c 


0.3116L 1 and 0.248814 6 


0.2489 ± 0.0002 


0.0198 ±0.0005 


0.282 ±0.005 


0.267 ± 0.005 


Pic 








0.0793 ± 0.0004 0.0992 ± 0.0005 


cr c 






0.3129 ±0.0005 






P^iPc) 




0.63 ±0.02 


0.34 ±0.02 


0.53 ±0.02 


0.54 ±0.02 


V 


0.875 


0.875 ± 0.070 


1.05 ±0.05 


0.875 ± 0.070 


0.875 ± 0.070 


P 


0.417 


0.43 ±0.04 


1.4 ±0.1 


0.43 ±0.04 


0.42 ±0.04 


T 






2.85 ±0.05 


2.19 ±0.05 





TABLE I: Comparison of the critical thresholds and critical exponents for conventional bond percolation, the loop percolation 
(LP) model and the loop Hamiltonian (LH) model. a refers to conventional site percolation and 6 refers to conventional bond 
percolation. 



constraint (jl]) of a Hamiltonian (||) . Our results are sum- 
marized in Tabic Q. 

We found that in 2d and 3d the LP model belongs to 
the universality class of the conventional (bond or site) 
percolation ||. A plausible explanation for this obser- 
vation is the following: We map an occupied (empty) 
plaquette onto a occupied (empty) site on an appropri- 
ate lattice. Fig. || illustrates the mapping in 2d: the two 
loops on a square lattice of system size L, Fig. || (left), 
are mapped onto two clusters consisting of 2 and 4 occu- 
pied sites on a (dual) square lattice of size L — 1, Fig. || 
(right). The resulting clusters of sites are clusters of fi- 
nite extended objects like k-mers in Ref. Q, which show 
the same universal behavior at the percolation transition 
as conventional site percolation fij. Also we found that 
the non-trivial orientation of the loops in the LP model is 
irrelevant for the universality class. We expect that this 
also holds for higher spatial dimensions d with the same 
critical dimension d c = 6 as for conventional percolation 
ill- 

For the loop Hamiltonian (LH) model Eq. O) we found 
evidence for an unconventional universality class of per- 
colation in 2d and 3d, the exponents are listed in table 
I. In 2d, we obtained a rather large correlation length 



exponent v = 3.3 ± 0.3 , which possibly indicates an 
infinite critical exponent for L — > oo known from the 
Kosterlitz-Touless (KT) phase transition fl6| . On the 
other hand, since our Hamiltonian (|^) has no XY, a 
Kosterlitz-Touless (KT) transition can be ruled out. In- 
deed, from applying a KT-form of finite-size scaling to 
our data, we could not find any acceptable data collapse. 
We like to remark that loop percolation with a large cor- 
relation exponent v also appears in integer Quantum Hall 
systems, where the loops represent equipotential lines 

& 

It would be interesting to study the universal behavior 
of the pure (i.e. a — 0) LH model (||) for finite tem- 
peratures T. In 3d, such a thermal-driven loop percola- 
tion phase transition has been studied Jj| for a different 
model. 
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